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ABSTRACT 

We report the detection of a weak X-ray point source coincident with the nucleus of the bulgeless disk 
galaxy NGC 3621, recently discovered by Spitzer to display high ionization mid-infrared lines typically 
associated with AGN. These Chandra observations provide confirmation for the presence of an AGN in 
this galaxy, adding to the growing evidence that black holes do form and grow in isolated bulgeless disk 
galaxies. Although the low signal-to-noise ratio of the X-ray spectrum prevents us from carrying out 
a detailed spectral analysis of the nuclear source, the X-ray results, combined with the IR and optical 
spectroscopic results, suggests that NGC 3621 harbors a heavily absorbed AGN, with a supermassive 
black hole (SMBH) of relatively small mass accreting at a high rate. Chandra also reveals the presence 
of two bright sources straddling the nucleus located almost symmetrically at 20" from the center. Both 
sources have 0.5-8 keV spectra that are well- fitted by an absorbed power-law model. Assuming they 
are at the distance of NGC 3621, these two sources have luminosities of the order of 10^^ erg s""'^, 
which make them ultra-luminous X-ray sources (ULXs) and suggest that they are black hole systems. 
Estimates of the black hole mass based on the X-ray spectral analysis and scaling laws of black hole 
systems suggest that the 2 bright sources might be intermediate mass black holes with A/bh of the order 
of a few thousand solar masses. However, higher quality X-ray data combined with multi-wavelength 
observations are necessary to confirm these conclusions. 

Subject headings: Galaxies: active - Galaxies: nuclei - X-rays: galaxies 



1. INTRODUCTION 

The well-known correlation between the black hole mass, 
Mbh, and the host galaxy stellar velocity dispersion cr^ 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000) implies 
that black hole growth and the build-up process of galaxy 
bulges are closely related. The fact that the vast ma- 
jority of AGNs reside in galaxies with prominent bulges 
(e.g. Heckman 1980; Ho, Fihppenko, & Sargent 1997; 
henceforth H97) suggests that bulges play a critical role 
in the formation and growth of nuclear black holes. How- 
ever, exceptions do exist, the most notable example being 
NGC 4395, a nearby Sd galaxy that lacks a bulge and 
shows optical and X-ray properties typical of a Seyfert 1 
galaxy (e.g., Fihppenko & Ho 2003; Moran et al. 2005) 
with an inferred black hole mass of 3.6 x 10^ Mq (Peter- 
son et al. 2005), much less massive than black holes found 
in galaxies with massive bulges. Until recently, this galaxy 
was an anomaly in that it was the only essentially bulgeless 
disk galaxy to harbor a nuclear black hole. It furthermore 
contains a nuclear star cluster (Fihppenko & Ho 2003), 
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providing a unique opportunity to study the connection 
between nuclear star clusters and black holes. More re- 
cently, Greene & Ho (2004; 2007) searched the Sloan Digi- 
tal Sky Survey for galaxies with similar intermediate mass 
black holes and found more examples of broad-line AGNs, 
many of which reside in disk-dominated galaxies that lack 
a classical bulge. Other more nearby examples of AGNs 
or AGN candidates in late-type galaxies have been found 
based on infrared, optical, and X-ray studies (Satyapal et 
al. 2008; Shields et al. 2008; Ghosh et al. 2008; Barth 
et al. 2008; Thornton et al. 2008; Dewangan et al. 2008; 
Desroches & Ho 2009). Although these studies demon- 
strate that AGNs reside in late-type disk galaxies, NGC 
4395 remains exceptional since it is both near enough for 
detailed study and is a truly bulgeless disk galaxy 

The investigation of AGNs in bulgeless galaxies is im- 
portant for the following reasons. First, it can shed light 
on the connection between black holes and galaxy forma- 
tion and evolution, by addressing the following outstand- 
ing questions: How does the black hole mass and accretion 
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rate relate to the properties of the parent galaxy in the case 
of bulgeless galaxies? Second, since bulgeless galaxies are 
thought to be in an early phase of their evolution, they al- 
low us to probe the mechanisms for growth of supermassive 
black holes which are still poorly known. Constraints on 
the low-mass end of the black hole mass function provide 
the best discriminant among models that predict the evo- 
lution of black holes from extremely high rcdshifts to the 
present time (e.g., Volonteri, Lodato, & Natarajan 2008). 
But before we can carry out any of the above tests, we 
have to locate the black holes themselves. 

The search for weak nuclear activity associated with 
black holes in extremely- late type normal galaxies is obser- 
vationally challenging because these galaxies are typically 
dusty and have significant star formation rates. Therefore, 
such a search is best carried out in energy bands that do 
not suffer from heavy dust or gas extinction and in which 
contaminating emission from star-forming regions is rela- 
tively weak. Thus, the mid-IR and hard X-ray bands are 
well suited for this purpose and, in fact, studies in these 
two bands complement each other well. 

Observations of the mid-IR fine structure lines provide 
us with one of the few definitive tools for discovering buried 
AGNs in dusty galaxies. AGNs show prominent high ex- 
citation fine structure lines whereas starburst and nor- 
mal galaxies are characterized by lower excitation spectra 
characteristic of Hii regions ionized by young stars (e.g. 
Genzel et al. 1998; Sturm et al. 2002; Satyapal et al. 
2004). One of the key prominent lines is the [Nev] 14/im 
line, which cannot be produced in Hii regions surround- 
ing young stars, the dominant energy source in starburst 
galaxies, since even hot massive stars emit very few pho- 
tons with energy sufficient for the production of this ion. 
The power of these diagnostics in finding buried AGNs 
has been convincingly demonstrated by the discovery of 
a population of AGNs in galaxies that display optically 
"normal" nuclear spectra (e.g., Satyapal et al. 2004). 

Once the presence of an AGN is ensured via IR spec- 
troscopy. X-ray observations represent the most effective 
means to investigate the AGN properties, since: 1) The 
X-rays are produced and reprocessed in the inner, hottest 
nuclear regions of the source. 2) The penetrating power of 
(hard) X-rays allows them to carry information from the 
central engine without being substantially affected by ab- 
sorption, and to probe the presence of the putative torus 
(provided that Nh < 10^* cm~^). 3) Compared to opti- 
cal and UV radiation. X-rays are far less affected by host 
galaxy contamination. 

Using high resolution IRS observations from the Spitzer 
archive, Satyapal et al. (2007, 2008) started to investi- 
gate the IR spectral properties of late-type (Sc or later) 
galaxies. Remarkably, of the 33 late- type galaxies with 
no definitive signatures of AGNs in the optical, 8 display 
[No v] emission. This result suggests that AGNs might be 
common not only in early-type galaxies but also in late- 
type ones, extending the ubiquity of the AGN phenomenon 
in the local Universe. In their sample of galaxies showing 
robust evidence for a [Nev] hne, NGC 3621 (z = 0.0024) is 
unique in that, like NGC 4395, it does not show evidence 
of any bulge component, it contains a nuclear star cluster 
(Barth et al. 2009) and is the closest one, with a distance 
of 6.6 Mpc bulgeless galaxies based on Cepheid measure- 



ments (Freedman et al. 2001). Very recently, Barth et 
al. (2009), using high quality optical spectra through a 
small aperture, detected a Seyfert 2 emission-line spec- 
trum, supporting the hypothesis that NGC 3621 harbors 
an accreting supermassive black hole. Additionally, the 
measurement of the stellar velocity dispersion provided an 
upper limit for the black hole mass of 3 x 10^ Mq. X-ray 
observations can strengthen the case for presence of an ac- 
creting black hole in this galaxy and allow us to determine 
some of the basic properties of the system, such as the ion- 
izing luminosity and the accretion rate, given an estimate 
of the black hole mass. 

Here, wc report the first X-ray study of NGC 3621 based 
on a 25 ks Chandra observation. The primary goal of this 
project is to verify whether this galaxy harbors an ac- 
creting supermassive black hole. As (explained below, our 
Chandra observation reveals the presc;nc;c! of 2 additional 
bright X-ray sources located ^20" from the center. There- 
fore, a secondary goal of this work is to investigate the 
nature of these 2 enigmatic X-ray sources. 

The outline of the paper is as follows. In § 2 we describe 
the observations and data reduction. The X-ray results are 
reported in § 3. In § 4 we discuss our findings, and in § 5 
we summarize the main results and draw our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 3621 was observed with Chandra ACIS-S on March 
6, 2008. The Chandra data, provided by the Chandra X- 
Ray Center, were processed using CALDB v. 3.4.2 following 
standard criteria. Only events with grades 0, 2-4, and 
6 and in the energy range 0.5-8 keV were retained. We 
also checked that no flaring background events occurred 
during the observations. The extraction of images, light 
curves, and spectra, as well as the construction of response 
matrices, was carried out following standard procedures. 
Background spectra and light curves were extracted from 
source-free regions on the same chip of the source. The 
counts were extracted from circular regions with radius of 
3". 

The spectral analysis was performed using the XSPEC 
V . 12 . 3 . software package (Arnaud 1996; Dorman fc Ar- 
naud 2001). For the brightest source, the spectra were 
re-binned in order to contain at least 15 counts per chan- 
nel, which is appropriate for the use of the statistic. 
For the second off-nuclear source, given the low num- 
ber of counts, we used the Cash statistic (Cash 1979) to 
fit the un- binned spectra over an energy range of 0.5- 8 
keV (in the observer's reference frame), where the cal- 
ibration is best known (Marshall et al. 2005) and the 
background negligible. The errors on spectral parameters 
are at 90% confidence level for one interesting parameter 
(AX^, AC = 2.71). 

3. RESULTS 

Chandra, with its sub-arcsecond spatial resolution, is the 
best X-ray satellite to investigate the spatial distribution 
of the X-ray emission. This is crucial for AGN studies 
and specifically for detection of weak AGNs, since these 
are typically associated with compact un-resolved regions 
as opposed to starburst activity, which instead is char- 
acterized by extended emission (e.g., Dudik et al. 2005; 
Gonzalez-Martm et al. 2006; Flohic et al. 2006). 
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Figure [T] shows the Chandra contours overlaid with the 
Spitzer IRAC 3.6/im image of the NGC 3621 galaxy from 
the SINGS survey (Kennicutt et al. 2003). The nuclear X- 
ray source shows a clear correspondence with the IR peak. 
Figure m shows the Chandra ACIS-S image of NGC 3621 
overlaid with 1.4 GHz contours obtained from a VLAQ 
snapshot observation at 20" resolution (160 s, archival 
program AC345). The nuclear X-ray source (hereafter 
source A, marked by a black box) is located at the center 
of the image. 

Interestingly, two bright off-nuclear X-ray sources are 
detected at ~20" from the center. While no clear IR as- 
sociation is evident for the 2 X-ray off- nuclear sources, 
the brighter X-ray source B may have a radio counter- 
part in the VLA map with a peak of 5.2 mJy/beam that 
is ^4" offset (to the north of) the X-ray source. At 
13=6.6 Mpc, the corresponding 1.4 GHz radio luminos- 
ity is 3 X 10^^ erg Hz^^ for an assumed point source. 
At the position of the X-ray nucleus and source C, the 
radio map is dominated by diffuse emission from the 
galaxy and the formal la upper limits are 3.9 mJy/beam 
and 4.9 mJy/beam (the radio surface brightnesses at the 
X-ray position), corresponding to point source limits of 
;S (2 - 3) X 10^'' erg s~^ Hz~^ A higher resolution VLA 
8.5 GHz snapshot (107 s, archival program AR555) did 
not reveal any significant emission at these locations with 
point source limits of 0.5 mJy (5cr) at 0.53" x 0.22" res- 
olution. The vastly differing resolutions of the two ra- 
dio maps preclude any practical constraints on the radio 
source spectra. 

At first glance the symmetric placement of the two off- 
nuclear sources relative to respect to the center of the 
galaxy, and their apparent association with radio emission, 
suggests that they may be related to the central source. 
They resemble hot-spots produced by the interaction of 
bipolar outflows with the interstellar medium. However, 
the results of our analysis of their time-averaged X-ray 
spectra (see below) seem to rule out this scenario. 

3.1. Source A: The Central Source 

The X-ray central region is located at 
RA=11'^18'"16.5P, DEC=-32° 48' 50.4" (the coordinates 
here and in the rest of the paper refer to the J2000 epoch, 
and have been determined using the celldetect task in 
CIAO). Note that the Chandra position of the central 
source is fully consistent with the NGC 3621 coordinates 
from the 2MASS catalog, RA=ll''18'"16.50", DEC=-32° 
48' 50.6" (Skrutskie et al. 2006). 

Using the CIAO tool dmextract, we extracted the X- 
ray counts from a central circular region of radius of 1.5" 
(black thick circle in Figure [3]) and from 8 background 
regions surrounding it, as shown in Figure [3l After back- 
ground subtraction, we detect 21 ± 6 X-ray counts from 
the source, which is significant at the 3.7 a level. The 
uncertainty was calculated using the approach of Gehrels 
(1986), which is appropriate for a small number of counts. 

Unfortunately, with such few counts, it is impossible 
to spectrally characterize source A; neither direct spectral 
fitting nor hardness ratio calculations yield useful informa- 
tion. Examining the morphology in several, narrow energy 

^ The National Radio Astronomy Observatory is operated by Assoc 
Science Foundation. 



bands only reveals that most of the counts (~70%) are in 
the energy range 0.8-1.8 keV, which coincides with the 
peak of the Chandra ACIS effective area curve. Therefore, 
we resort to deriving some estimates of the absorbed flux 
and corresponding intrinsic luminosity using the simula- 
tion tool PIMMS. 

Assuming a power-law model with photon index F rang- 
ing from 1.7 to 2.0 (which is typical for AGNs) and a 
Galactic column density of A^h = 5.74 x lO^^cm"^, we ob- 
tain an absorbed flux of -Fo.s-skcV ^ 6.5 x 10^^^ erg cm~^ 
s-i (F2_iokcV 4 X 10-1^ erg cm'^ s'^). The Seyfert 
2 spectroscopic classification and the relatively large bolo- 
metric luminosity, inferred by both IR and optical studies 
(Satyapal et al. 2007; Barth et al. 2009), suggest the pres- 
ence additional absorption, local to NGC 3621. In order 
to infer plausible values for the intrinsic luminosity, we 
assumed some additional intrinsic absorption with column 
densities ranging from 5 x 10^^ to 5 x 10^^ cm~^. The lower 
value of A^H is consistent with the one derived from spec- 
tral fitting of the 2 off- nuclear sources (see below), whereas 
the upper value is a typical value found in non Compton- 
thick Seyfert 2 galaxies (see, e.g., Figure 3 of Bassani et 
al. 1999). The resulting values of the estimated luminosi- 
ties, which mostly depend on the choice of intrinsic 
(variations of F cause luminosity differences of the order 
of a few percent only) are Lo.s-skcV ~ 7 x 10'^'' erg 
(L2-iokoV ~ 5 X 10^^ erg s^^) and Lo.5-8keV ~ 3 x 10^^ 
erg s~^ (i2-iokcV ^ 2 X 10'^^ erg s~^), for the moderate 
and heavy absorption scenarios, respectively. 

In summary, the Chandra observations reveal the pres- 
ence of a weak point-like X-ray source at the center of 
NGC 3621, which confirms the presence of an AGN. How- 
ever, the paucity of counts hampers the study of the nature 
of this source and more specifically whether this system is 
heavily absorbed or not. 

3.2. Source B: the Brightest Off- Nuclear Source 

Source B is fairly X-ray bright (more than 2000 counts 
were collected during the 20.6 ks net exposure) located at 
RA=lll^l8'"15.16^ DEC=-32° 48' 40.6". Since Source B 
is much brighter than the other sources in the Chandra 
field of view, one may have the impression that its emis- 
sion is spatially extended. However, a comparison of the 
brightness radial profile with the ACIS point spread func- 
tion reveals that X-ray spatial distribution is fully consis- 
tent with that of a point-like source. 

Figure H shows the Hubble Space Telescope {HST) ACS 
WFC images centered around the position of Source B. 
The position of the off-nuclear X-ray source was observed 
at the edge of the field of HST images obtained in the 
F435W, F555W (not shown), and F814W filters (each with 
3x360s exposures; program 9492) on 2003 February 16. 
Aligning the X-ray and optical nucleus in the Chandra 
and i/ST images, we find a faint optical source 0.33" from 
the center of the X-ray position of source B, although only 
the eastern ~l/2 of the X-ray error circle was covered in 
the WFC images. Ulitizing r=0.2" circles, we measured 
count rates for the optical source including four surround- 
ing source-free regions for background estimation. The 
optical source is a 4.5(t detection in the F435W image, a 
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lower significance detection (2.8(7) in the F555W image, 
and is undetected in the F814W image, where 1 a errors 
were determined by calculating the standard deviation of 
the background count rates. Count rates were converted to 
flux density utihzing the PHOTFLAM and PHOTPLAM 
keywords in the FITS header and the aperture corrections 
tabulated in Sirianni et al. (2005), and are reported in 
Table 1. 

The relatively high X-ray count rate of Source B allows 
us to investigate not only the X-ray spectral properties 
via direct fitting of the spectrum but also the temporal 
properties. We have extracted the light curve in the 0.5-8 
keV energy range, from a circular region with a radius of 
3" centered on source B. This light curve is shown in Fig- 
ure [HI Although there appear to be some low-amplitude 
fluctuations toward the end of the observation, according 
to a test the light curve is consistent with the hypoth- 
esis that the source flux is constant during the Chandra 
pointing (xVdof=24.20/22, P^2 = 0.34). 

The X-ray spectrum of source B, obtained from the 
same extraction region as the light curve, is poorly fit- 
ted by a simple power law with Galactic absorption 
(x^/dof=154. 12/112), showing a deficit at low energies in- 
dicative of additional absorption. Indeed, by leaving the 
column density of the absorber free to vary, a much better 
fit is obtained (x^/dof=50. 47/111) with the following pa- 
rameters: iVn = (3.4 ±0.6) X lO^i cm^^ and F = 2.7 ±0.3. 
The spectrum of source B, with the best fit overlaid, 
and the data-to-model ratio are shown in Figure [6l The 
corresponding absorbed flux is -Fb.5-8kcV = 7 x 10^ ^'^ 
erg cm~^ s~^ and the intrinsic luminosity, assuming that 
source B is located at the same distance as NGC 3621, 
-^^o.5-8kov = 7 X 10^^ erg s'^ (La-iokeV = 2 x 10^^ erg 
s"-'^). We have also tried to add a thermal component (de- 
scribed by the XSPEC model apec). Although the over- 
all fit remains acceptable (x^/dof=48. 97/108), the spec- 
tral parameters of the thermal component {kT and abun- 
dances) remain largely unconstrained. 

The large luminosity, well in excess to the Eddington 
limit for a neutron star, suggests that source B can be a 
black hole system. As a consequence, the X-ray emission 
is likely to be produced by the Comptonization of seed 
photons emitted by the underlying accretion disk. It is 
therefore instructive (and indeed a necessary step to con- 
strain Mbh in this system with the method described in 
§4.2) to fit the Chandra spectrum of source B with the 
physically motivated BMC model instead of a phenomeno- 
logical power law (see Titarchuk et al. 1997 for further 
details). 

The Bulk Motion Comptonization model (BMC) is a 
generic Comptonization model able to describe equally 
well bulk motion and thermal Comptonization, although 
it was historically developed to describe the Comptoniza- 
tion of thermal seed photons by a relativistic converging 
fiow (Titarchuk et al. 1997). The BMC model is char- 
acterized by 4 free parameters: the temperature of the 
thermal seed photons kT, the energy spectral index a, 
a parameter \og{A) related to the Comptonization frac- 
tion / by the relation / = A/(1 + A) (where / is the ra- 
tio of Compton scattered photons over the seed photons), 
and the normalization A^bmc- In simple words, the BMC 
model convolves the thermal seed photons and a generic 



Comptonization Green's function producing a power law. 
As a consequence, this model adequately fits X-ray spec- 
tra of accreting black holes, which at the zeroth order are 
always characterized by a hard power law, which is widely 
believed to be produced by the Comptonization of a seed 
thermal component. In addition to being a comprehensive 
Comptonization model, the BMC model presents 2 impor- 
tant advantages with respect to the power law model (PL): 

1) Unlike the PL, which is a phenomenological model, the 
BMC parameters are computed in a self-consistent way; 

2) Unlike the PL, the power law produced by BMC does 
not extend to arbitrarily low energies and thus does not 
affect the normalization of the thermal component nor the 
amount of local absorption, which is often present around 
accreting objects. 

The resuhing best fit (xVdof=49.90/110), which is sta- 
tistically indistinguishable from the PL fit, yields the 
following parameters: A^h — (3.1 ± 0.1) x 10^^ cm~^, 
kT = 0.1 ± 0.1 keV (the temperature of the seed pho- 
tons), a — 1.6 ± 0.3 (the spectral index, which is related 
to the photon index by F = a + 1), log A — 2 (fixed at the 
best fit value), and A'^bmc = l-2lo 4 ^ 10"^, which is in 
units of (L/10^^ erg s~^)(10 kpc/d)^ with L and d being 
the luminosity and distance of the object, respectively. 

3.3. Source C: the second off-nuclear source 

The measured X-ray position of Source C is 
RA=llhl8■"18.23^ DEC=-32° 48' 53.0". Figure [7] shows 
the HST ACS WFC images (described in the previous sub- 
section) centered around the position of Source C. A faint 
optical source is detected 0.37" away from the Chandra 
position in all 3 filters. An additional redder source is 
visible 0.5" away but detected only in the F814W filter 
image. Additional ACS images from the HST same pro- 
gram obtained on Feb 3 confirm the two optical sources. 
The photometric results from this latter dataset are within 
the la errors with the exception of the F435W measure- 
ment of CI where we measure a flux that is 2a lower than 
that reported in Table 1. 

The moderately low X-ray count rate (~200 counts col- 
lected during the net 20.6 ks exposure) allows a direct 
fltting of the un-binned spectrum using the C— statistics, 
which can be used to estimate parameters values and con- 
fidence regions but does not provide a goodness-of-fit. 
To this aim, on can use the XSPEC command goodness 
that performs Monte Carlo simulations of spectra based 
on the chosen model. This procedure yields the per- 
centage of simulations with the fit statistic lower than 
that for the original data. If the best-fitting model is a 
good representation of the data the percentage should be 
around 50%; values ^50% indicate that the data are over- 
parameterized, whereas values close to 100% indicate that 
the fit is poor. 

Source C appears to be reasonably well fitted with 
a simple power-law model with absorption in excess to 
the Galactic A'h (see Figure [5]), although the goodness 
command yields 74% indicating that this model is only 
marginally acceptable. The best fit parameters are: A^h = 
(5.ltij)xl02i cm-2 andF = 1.6±0.4. The corresponding 
absorbed fiux is Fo.5_8koV = 1-1 x 10^^'^ erg cm~^ s~^ and 
the intrinsic luminosity, assuming that source C is located 
at the same distance as NGC 3621, Lo.s-skoV = 7 x 10'^* 
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erg s 1 (L2_iokeV = 6 x 10^^ erg s 

Also in this case we have tried the Comptonization BMC 
model instead of a phcnomenological power law. The re- 
sulting best fit parameters are: A^h = (5.0 ± 2.2) x 10-^^ 
cm-2, kT = Omto^l keV, a = 0.5 ± 0.2 , log A = -1.15, 
A'^BMC = (6.3 ± 5.0) X 10~®. The goodness command 
(98%) indicates that this model is statistically worse than 
the power law. However, as explained above, the results 
from the fitting of the BMC model represent the first neces- 
sary step to attempt an estimate of Mbh in this putative 
black hole system. 

4. DISCUSSION 
4.1. The central source 

The primary goal of our X-ray investigation is to ver- 
ify whether NGC 3621 hosts a buried AGN, as suggested 
by recent IR and optical spectroscopic studies. The pres- 
ence of a weak X-ray source coinciding with the center 
of the galaxy seems to be consistent with this scenario. 
Unfortunately, however, the very low count rate severely 
hampers our investigation, preventing any X-ray spectral 
characterization of this source. Nevertheless, combining 
the Chandra results with findings from recent IR and op- 
tical spectroscopic studies, wo may address several impor- 
tant questions. 

Before proceeding further, it is worth summarizing the 
main results from longer wavelength studies that will be 
used in our discussion. From the IR study, we will use 
the bolometric luminosity, Lboi = 5 x 10^^ erg s~^, which 
was derived from the i[NeV] ~ -^'boi correlation (Satyapal 
et al. 2007). From the optical investigation of Barth et 
al. (2009) we will make use of the following information: 
1) NGC 3621 is spectrally classified as a Seyfert 2; 2) The 
measured [O iii] A5007 flux is f[oiii] = (22.4±0.4) x lO'^'^ 
org cm~^ s~^, which, based on inspection of the two- 
dimensional optical lino ratios, probably encompasses only 
10% of the narrow line region (NRL). Note that if the last 
hypothesis is correct, the bolometric luminosity derived by 
■^[O III] ~ -^boi correlation is fully consistent with the value 
inferred from Z/[N(,v]; if not, iboi derived from L[oiii] is 
lower by a factor of ^10. 

The first outstanding question is whether the central 
source is heavily absorbed or intrinsically weak. Based on 
the Chandra, coimt rate, the 2 10 kcV luminosity we have 
inferred has a range of over 2 orders of magnitude (from 
a few xlO"*^ to a few xlO'^^ erg s^^), depending on the 
adopted intrinsic A^h (5 x 10^^ and 5 x lO^'^cm"^, respec- 
tively). 

We can estimate the -L2-10 koV luminosity using the 
LiNeV] luminosity assuming the L[^ey]/L2~io koV ratio 
in standard AGN. Using the limited number of standard 
AGN with both [NeV] and 2-10 keV observations (based 
on recent Spitzer observations of AGN from Dudik et al. 
2007 and Gorjian et al. 2008), we find an average ratio of 
L[NeV]/ L2-10 keV — 10^"^. This result, combined with the 
estimated value of L[NeV] for NGC 3621, ~ 5 x 10^'^ erg 
s^^, predicts a 2 10 keV luminosity of the order of lO''" 
erg s^^, which is more than two orders magnitude greater 
than the observed luminosity in the weakly absorbed sce- 
nario. Therefore, this result strongly suggests that NGC 
3621 is heavily absorbed. 



A second indirect way to test the value of Lx employs 
the correlation between Lx — -L[oiii]- Using the results 
from Mulchaey et al. (1994) for a sample of Seyfert galax- 
ies, and, more specifically, utilizing the best-fit linear cor- 
relation derived by their Figure 3c, we obtain a luminosity 
of Lx = 9.8 x 10^^ erg s~^, which is in broad agreement 
with the heavily absorbed scenario. 

A third alternative way to constrain Ajj (and hence Lx) 
is based on the -Fx/-F[oiii] — A^h diagram proposed by 
Guainazzi et al. (2005) and derived from the work of Bas- 
sani et al. (1999) on Seyfert 2 galaxies. In this diagram, 
which shows a clear anti-correlation between fx/L'[oiii] 
and A^H, -Fx is the absorbed 2-10 keV flux whereas -^[0111] 
is corrected for absorption based on the Ha/H(3 ratio (see 
the appendix of Bassani et al. 1999 for details). Depending 
on whether -/^[oiii] encompasses 10% or the entire [O iii] lu- 
minosity of the NLR, the Fx/-F'[oiii] ratio lies between 0.1 
and 1, which in turn translates into A^h ranging between 
~ 10^'^ and ^ 10^** cm^^. This result again lends support 
to the heavily absorbed scenario. Finally, this scenario 
is supported by the direct comparison between Lx and 
Lboi- Indeed, if we use the value obtained with moderate 
absorption, 5 x 10^''erg s~^, we derive a bolometric correc- 
tion factor of K2-iokcV = Lboi/L2-iokcV > 1000, which 
is unreasonably high and inconsistent with any reasonable 
spectral energy distribution (SED) for AGNs. 

Adopting the heavily obscured scenario, we can now try 
to constrain the most important parameters for a black 
hole system, namely the accretion rate, m, and the black 
hole mass, Mbh- Assuming that NGC 3621 has an X-ray 
luminosity of few xlO^^ erg s~^, the bolometric correction 
factor becomes K;2-iokeV ^100. This result directly trans- 
lates into an estimate the Eddington ratio value (and hence 
on m) by virtue of the existing positive correlation between 
the bolometric correction factor K;2-iokeV and -Lboi/-£'Edd 
(Vasudevan and Fabian 2009). According to Vasudevan 
and Fabian (2009), this bolometric correction is character- 
istic of systems accreting at higher rate: LboiZ-^^Edd ^ 0.2. 
This finding in turn can be used to constrain the black 
hole mass: Assuming that Lboi = 5 x 10^^ erg s~^and 
Lho\/L-Edd 0.2, we obtain Mbh ^ 2 x 10^ M©, which is 
a fairly low value and apparently below the extrapolation 
of the Mbh — cr correlation (see, e.g., Fig. 10 of Barth et al. 
2009). The small Mbh coupled with the high rn appears to 
be in line with the scenario where bulgeless galaxies host 
small supermassive black holes in their early stage that 
undergo a phase of very vigorous accretion, and is consis- 
tent with recent findings from Greene & Ho (2007) and 
Barth et al. (2008) on Seyfert 1 and Seyfert 2 galaxies, 
respectively. 

We should point out that the scenario of a heavily ab- 
sorbed AGN in NGC 3621 is not without caveats. The 
correlation between the [NeV] luminosity and the bolomet- 
ric luminosity in standard AGN shows a scatter 0.44 dex, 
which corresponds to a factor of ^3 uncertainty in Lboi 
derived by L[ncV] (Satyapal et al. 2007). Furthermore, 
its applicability at this luminosity range is uncertain. In 
addition the fraction of the [QUI] luminosity encompassed 
by the optical spectrum obtained by Barth et al. (2009) 
and attributable to the AGN is also uncertain. However, 
taken collectively, there is strong support for the hypothe- 
sis that NGC 3621 harbors a buried AGN with luminosity 
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of the order of a few times 10^^ erg s ^. 

4.2. The off-nuclear sources 

The secondary goal of this work is to investigate the na- 
ture of the 2 bright off-nuclear sources detected by Chan- 
dra. Based on the X-ray luminosity values derived from 
the spectral analysis {L2-iokeV = 2 x 10'^^ and 6 x 10"^* 
erg s~^ for sources B and C, respectively), with the as- 
sumption that these sources are located within NGC 3621 
(i.e., at a distance of 6.6 Mpc), we can already derive the 
first important conclusion: these are ultraluminous X-ray 
sources and their compact objects are likely black holes, 
since their X-ray luminosity exceeds the Eddington limit 
for neutron stars. 

The lack of strong optical/UV counterparts (see Table 
1), the fact that they have similar intrinsic absorption that 
is well in excess with respect to the Galactic value, and the 
apparent association with the NGC 3621 radio emission 
suggest that these sources are not located at significantly 
smaller distances than NGC 3621. 

However, a priori, we cannot rule out that they are back- 
ground AGN that appear to be located close to the cen- 
ter of NGC 3621 just by chance. To test this hypothesis, 
we have performed the following test. Utilizing the re- 
sults of the Chandra spectral analysis we have derived the 
monochromatic X-ray luminosity I2 kcV for different val- 
ues of redshift ranging between z = 0.05 — 1. From the 
well-known tight correlation I2 krv — I ° for standard 

° 2500A 

AGN (we used equation Ic from Steffen et al. 2006) we 
have determined the corresponding values of I ° and 

2500A 

the relative fluxes. Finally, the expected values for the 
flux at the wavelengths probed by the HST observations 
have been obtained by converting the 2500A fluxes assum- 
ing a typical slope of 0.7 {fi, oc i/"" ''). 

For Source B the predicted optical flux ranges between 
230-370 fiJy (depending on the optical band) at 2; = 0.05 
and 6-9 mJy at z = 1. For the weaker X-ray Source C, we 
obtain optical fluxes of the order 14-21 /iJy at z = 0.05 
and 345-535 /iJy at z = 1. Since these expected values for 
standard AGN are at least 2 orders of magnitude larger 
than the measured optical fluxes reported in Table 1, we 
can rule out the hypothesis of background AGN. We also 
applied extinction corrections by assuming the absorption 
column density derived from the X-ray spectral fits, A^h, 
and converting it into optical extinction with the relation 
A{V)/Nn = 5.3 X 10-22 (Cox & Allen 2000). This cor- 
responds to A{V) —1.8 and 2.7, which in turn translates 
into flux correction factors of ~5 and '^10 for source B and 
C, respectively. Nevertheless, the measured optical fluxes 
are still inconsistent with the predicted values for standard 
AGN. 

We can therefore reasonably assume that both off- 
nuclear sources reside in NGC 3621. At a distance of 6.6 
Mpc, the measured fluxes of their putative optical coun- 
terparts correspond to optical luminosities of the order 
of ^ 10'^^ — 10'^^ erg s"^, which, combined with the esti- 
mated X-ray luminosities, appear to be in line with typical 
values of Lx/-^opt measured in X-ray binaries (e.g., Rit- 
ter & Kolb 2003). In the hypothesis that these sources 
are indeed hosted by NGC 3621, it is important to deter- 
mine whether they are analogs of Galactic microquasars 
like GRS 1915-^105, or intermediate black holes (IMBHs) 



close to the center of the galaxy. 

In order to shed light on the nature of these sources, and 
specifically derive useful constraints on their Mbh and dis- 
tance, we will make use of two different methods, that are 
based on the X-ray spectral fitting results with the BMC 
model. Both methods rely on the generally accepted as- 
sumption that the X-ray spectrum is produced by Comp- 
tonization of a thermal component associated with an ac- 
cretion disk. The first method assumes a specific model for 
the accretion disk based on the seminal work from Shakura 
& Sunyaev (1973) and interpret kT as the color tempera- 
ture of the accretion disk (which is related to the effective 
disk temperature by a "hardening factor" Th). On the 
other hand, the second method only assumes that the un- 
derlying physics in black hole systems is the same irrespec- 
tive of the scale. Therefore, with the second method Mbh 
can be directly derived by scaling self-similar relationships 
between spectral (and temporal) parameters. For further 
details, see Shaposhnikov & Titarchuk (2009). 

The first method has been introduced by Shrader & 
Titarchuk (1999) and successfully applied to derive Mbh 
in several Galactic black holes (GBHs), in intermediate 
BHs, and also in a sample of Narrow Line Seyfert 1 galax- 
ies (see Shrader & Titarchuk 2003 for details). The sec- 
ond method, on the other hand, is based on the extension 
of the scaling technique presented by Titarchuk & Fiorito 
(2004) and Shaposhnikov & Titarchuk (2007). Specifically, 
it is based on the universal scalable relationship between 
the photon index and the normalization of the BMC model, 
which appears to be characteristic of all GBHs during their 
spectral transitions (Shaposhnikov & Titarchuk 2009). 

In simple terms, the second method can be summarized 
by the following steps. 1) Construct a F — A'^bmc pfot for 
a GBH of known mass and distance, which will be used as 
reference (hereafter denoted by the subscript r). 2) Com- 
pute the normalization ratio between the target of interest 
(hereafter denoted by the subscript f) and the reference ob- 
ject A^BMC,t/A'BMC,r at the corresponding value of F. 3) 
Derive the black hole mass using the following equation 

MBH,t = Mbh,!- X (A^BMC,t/A^BMC,r) X {dtldrf X /g (1) 

where MBH.r is the black hole mass of the reference object, 
A^BMC,t and A^BMC,r are the respective normalizations for 
target and reference objects, dt and d^ are the correspond- 
ing distances, and fa = cos 9r/ cos 6*4 is a geometrical fac- 
tor that depends on the respective inclination angles. The 
above formula is readily obtained by considering that a) 
the normalization is a function of luminosity and distance: 
A'bmc oc L/d^; b) the luminosity can be expressed by 
L (X t^Mbhw, where rj is the radiative efficiency, and c) 
assuming that different sources in the same spectral state 
(defined by the photon index) have the same 77 and to. 

To illustrate this method, in Figure [9] we show the 
r — A^BMC diagram for GRO J1655-40, which is the pri- 
mary reference source since the parameters of this system 
are the most tightly constrained: Mbh/Mq = 6.3 ± 0.3, 
i = 70° ±1°, d = 3.2±0.2kpc (Greene et al. 2001; Hjellm- 
ing & Rupen 1995; but see also Foellmi et al. 2006 for a 
discording view on the distance of GRO J1655-40). 

Since, unlike the first method from Shrader & Titarchuk, 
this method is relatively new and has never been tested 
before for AGNs or IMBHs, we apply it first to an AGN 
of known mass, as a check. We chose the broad-line 
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radio galaxy (BLRG) 3C 390.3 (z=0.056 corresponding 
to a luminosity distance of 247 Mpc), because, it has a 
well constrained mass as well as a well constrained disk 
inclination. One black hole mass estimate of Mbh = 
(5 ± 1) X 10^ M0 is based on the velocity dispersion (Nel- 
son et al. 2004; Lewis & Eracleous 2006) and another 
estimate of Mbh = (2.9 ± 0.6) x 10^ Mq based on re- 
verberation mapping (Peterson et al. 2004)- An incli- 
nation angle of 25° < i < 35° is derived from the radio 
data (Eracleous, Halpern & Livio 1996, Giovannini et al. 
2001, Lewis et al. 2005), while an angle of 26°t2o is de- 
rived from the profiles of the double-peaked optical emis- 
sion lines (Eracleous & Halpern 1994). Fitting the XMM- 
Newton EPIC pn spectrum of 3C 390.3 with the BMC model 
(and a Gaussian line to account for the Fe Ka emission), 
we obtain: kT = 0.106 ± 0.003 keV, a = 0.75 ± 0.01 , 
log(A) = 0.66 ± 0.03, and TVbmc = (8.7 ± 0.2) x 10"^. 
From Figure [H we infer that at F — 1.75 the reference 
normalization is A'^BMC.r — 0.3. Inserting this result and 
the corresponding quantities in equation (1) we obtain 
Mbh ^ 1-5 x 10** Mg. Using dGROJi655-4o = 1-7 kpc, 
as suggested by Foellmi et al. (2006), the resulting mass 
increases by a factor ^--^3.5, yielding Mbh — 5 x 10® Mq, 
which is in good agreement with the accepted value. 

Since the value of Mbh derived for 3C 390.3 is reason- 
ably close to the estimates obtained from optical data, 
we can try to apply this procedure (hereafter method 2 
for simplicity) and the Shrader & Titarchuk method 
(method 1) to the 2 X-ray enigmatic sources detected by 
Chandra. Using the results from the spectral fitting with 
the BMC model, for source B we obtain Mbh — (1-4 x 10^ - 
1.5 X 10*) M0 and Mbh ^ x 10^ - 3.4 x 10^) Mg, us- 
ing method 1 and method 2, respectively. The ranges of 
Mbh reflect the relatively large uncertainties of the spec- 
tral parameters. Similarly, for source C, the two methods 
yield respectively Mbh ^ (1-8 x 10^ - 1.3 x lO"*) Mg and 
Mbh ^ (2.3 x 10^ - 2 x 10'') Mg. 

These findings are important for a number of reasons. 
First, the consistency between the BH mass estimates ob- 
tained with 2 different methods lends support to the va- 
lidity of the 2 methods. Second, since the methods are 
independent, the mass obtained from one method can be 
used in the second method to verify the distance of the 
object, and the similarity between A/bh seems to confirm 
the hypothesis that the 2 enigmatic sources are located 
in NGC 3621. Finally, and perhaps most importantly, 
our analysis suggest the presence of 2 IMBHs within a 
few hundred parsecs from the center of the galaxy, which 
hosts a small SMBH. This result, if confirmed, may have 
important implications for our understanding of the forma- 
tion process of supermassive black holes. Indeed, several 
theoretical models predict the existence of off-nuclear su- 
permassive black holes, which eventually merge with the 
central SMBH (see, e.g., Volonteri & Madau 2008 and ref- 
erence therein). 

5. CONCLUSIONS 

We have analyzed a 25 ks Chandra observation of the 
bulgeless galaxy NGC 3621 to investigate the basic X-ray 
properties of the putative AGN and of 2 bright X-ray 
sources detected ~20" away from the center. The main 
results can be summarized as follows: 



• A weak X-ray source was detected at the center of 
the galaxy, lending support to the finding from IR 
and optical spectroscopic studies that NGC 3621 
hosts an AGN. Unfortunately, the signal-to-noise 
ratio hampers any spectral study of the source. We 
therefore used PIMMS to infer the possible lumi- 
nosity associated with the putative AGN. Assum- 
ing an absorbed power-law model with F = 1.7 — 2 
and TVh = 5 x 10^^ — 5 x 10^'^ cm~^, the resulting 
luminosity is io.5-8kcV ^ 7 x 10^'' — 3 x 10^^ erg s"-' 
(1/2-10 kov ~ 5 X 10^''-2 X 10^^ erg s^^), depending 
on the choice of intrinsic absorption. 

• Combining the X-ray luminosity with information 
derived from recent IR and optical observations, 
and using several independent tests, we infer that 
a heavily absorbed AGN scenario is more likely for 
NGC 3621. If this is the case, exploiting the di- 
rect correlation between bolometric correction fac- 
tor and Eddington ratio, we speculate that the 
AGN is in its early phase with a supermassive black 
hole of relatively small mass (A/bh — 2 x 10*^ Mg) 
accreting at high rate (iboi/iEdd ^ 0.2). 

• The Chandra observation also reveals the presence 
of 2 bright off-nuclear sources located ~20" away 
from the central source. Direct fitting of the X-ray 
spectra indicates that both are well described by 
absorbed power laws with luminosities of the or- 
der of 10'^^ erg s~^, in the hypothesis that they are 
located in NGC 3621. 

• The lack of strong optical/UV counterparts and the 
intrinsic absorption argue against the 2 off-nuclear 
sources being Galactic objects that appear to be 
in NGC 3621 by chance superposition. Similarly, 
the combination of X-ray and optical measurements 
seems to rule out the hypothesis of background 
AGN. Assuming that these sources are at the dis- 
tance of NGC 3621, leads to the conclusion that 
their luminosities exceed the Eddington limit for a 
neutron star and suggests that these are black holes. 

• In order to discriminate between 2 competing 
scenarios for these putative black hole systems 
(namely, IMBHs, or bright microquasars similar to 
GRS 1915-hl05 in NGC 3621), we have applied 
2 different methods that allow to determine the 
MBH-to-distance ratio, once the spectrum is fitted 
with the BMC model Comptonization model. Both 
methods consistently suggest that the 2 off-nuclear 
sources are located in NGC 3621 and have Mbh of 
a few thousand solar masses. 

In conclusion, our Chandra observation has confirmed 
that NGC 3621 is an extremely interesting object. 
Not only do the X-ray data support the scenario that 
NGC 3621 hosts a buried AGN, but they also suggest the 
possible presence of 2 IMBHs located relatively close to 
the center of the galaxy. These findings may play a crucial 
role in our understanding of the connection between galaxy 
formation and SMBHs at their center, and of the process 
leading to the formation of the SMBH itself. It must be, 
however, kept in mind that these results were derived from 
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data with limited statistics. Therefore, deeper X-ray expo- 
sures from more sensitive instruments complemented with 
observations at longer wavelengths (with emphasis on the 
radio band that is largely unaffected by absorption) will 
be necessary to confirm and possibly further our main con- 
clusions. 
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Table 1 

Possible optical counterparts of off-nuclear X-ray sources 



Source 


ra 


DEC 


Offset 


F„ (F435W) 


(F555W) 


F„ (F814W) 




(J2000) 


(J2000) 


(") 




(p3y) 


(MJy) 


B 


11:18:15.180 


-32:48:40.39 


0.33 


1.8 (2.5) ± 0.4 


1.4 (1.8) ± 0.5 


<4.5 (<5.2) 


CI 


11:18:18.248 


-32:48:53.29 


0.37 


0.52 (0.72) ± 0.06 


0.46 (0.59) ± 0.06 


0.74 (0.85) ± 0.16 


C2 


11:18:18.252 


-32:48:52.56 


0.52 


<0.19 (<0.26) 


<0.17 (< 0.22) 


0.52 (0.60) ± 0.16 



The limits of the monochromatic fluxes are 3 sigma. The pivot wavelengths of the filters are respectively 4318 A (for F435W), 5360 A (for 
F555W), and 8060 A (for F814W). Extinction corrected (based on Schlegel et al. 1998 maps) values are indicated in parentheses. 
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Fig. 2. — Chandra ACIS-S image of NGC 3621 in the 0.5—8 keV band with 1.4 GHz contours from a short VLA observation overlaid. The 
center of the galaxy, source A, is shown by the black box. Two off-nuclear bright sources, source B and source C, are located at ~20" from 
the center. 
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Fig. 3. — Chandra ACIS-S image of NGC 3621. The extraction region for the central source is indicated by the black thick circle. The 
surrounding circles are the background regions. 




Fig. 4. — H3T ACS WFC images (in 2 of the 3 observed bands) centered around the position of Source B, the brightest X-ray off-nuclear 
that is indicated with r=0.5" dashed circles. The images were smoothed with a Gaussian with a = 2 pixels (0.098"). The edge of the WFC 
field of view is obvious in all images with only about ~l/2 of the Chandra source B circle covered. Detected optical sources tabulated in 
Table 1 are indicated with solid r=0.2" circles. 




5000 10^ LSxlO-" 2x10* 

Time (s) 

Fig. 5. — Chandra ACIS-S light curve of source B in the 0.5-8 keV energy band; time bins are 1000 s. 



A Chandra view of NGC 3621 



11 



0.01 




Energy (keV) 

Fig. 6. — Spectrum of source B and data/model ratio to a simple power-law model modified by intrinsic photoelectric absorption. 
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Fig. 7. — i?5T ACS WFC images (in 2 of the 3 observed bands) centered around the position of Source C, indicated with r=0.5" dashed 
circles. The images were smoothed with a Gaussian with a = 2 pixels (0.098"). Detected optical sources tabulated in Table 1 are indicated 
with solid r=0.2" circles. 
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Fig. 8. — Spectrum of source C and data/model ratio to a simple power-law model modified by intrinsic photoelectric absorption. 
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Fig. 9. — Photon index plotted versus the normaUzation of the BMC model for the microquasar GRO J1655-40 during the decay of the 
2005 outburst (adapted from Shaposhnikov & Titarchuk (2008). The thick continuous line represents the best fitting function of the spectral 
trend. The dashed lines illustrate how we determine the value of A^BMCr that we use in equation (1) to estimate the black hole mass of the 
BLRG of 3C 390.3. 



